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Abstract 

The low-frequency magneto-optical properties of bilayer Bernal graphene are stud- 
ied by the tight-binding model with four most important interlayer interactions taken 
into account. Since the main features of the wave functions are well depicted, the 
Landau levels can be divided into two groups based on the characteristics of the 
wave functions. These Landau levels lead to four categories of absorption peaks in 
the optical absorption spectra. Such absorption peaks own complex optical selection 
rules and these rules can be reasonably explained by the characteristics of the wave 
functions. In addition, twin-peak structures, regular frequency-dependent absorption 
rates and complex field-dependent frequencies are also obtained in this work. The 
main features of the absorption peaks are very different from those in monolayer 
graphene and have their origin in the interlayer interactions. 
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Ever since few-layer graphenes (FLGs) have first been produced by mechanical exfoliation 1 ' 2 
and thermal decomposition, 3 ' 4 they have attracted considerable experimental and theo- 
retical research because of their exotic physical properties resulting from the hexagonal 
structure. Monolayer graphene (MG) and bilayer Bernal graphene (BBG; AB-stacked bi- 
layer graphene), two of the simplest few-layer graphenes, display very different electronic 
structures. MG exhibits isotropic linear bands near the Fermi level E F = 0; these bands 
become gradually anisotropic parabolic bands as the energy exceeds the region of ±0.5 
eV. 5 In BBG, the linear bands change into two pairs of parabolic bands due to the present 
interlayer interactions and stacking configurations. 6-14 Each pair is comprised of the oc- 
cupied valence and unoccupied conduction bands, which are asymmetric about Ep. 6 ' 9 The 
conduction and valence bands of the first pair slightly overlap about E F , while those of 
the second pair rise at the energies 12 0.34 eV (~ 71 + 76 ) and -0.39 eV (~ —71 + 76), 
respectively, where ^'s 8 ' 10 ' 15-17 (s—1, 3, 4 and 6) are four important interlayer interactions 
in BBG. The unusual electronic structures of MG and BBG have been experimentally ver- 
ified through angular- resolved photoemission spectroscopy. 18 ' 19 The main features of the 
zero-field energy bands, e.g., the broken linear dispersions in MG and the two groups of 
asymmetric electronic states in BBG, should be reflected in the features of the Landau 
levels (LLs) resulting from the magnetic fields. 12 

In the presence of a uniform perpendicular magnetic field B = B z, the zero-field en- 
ergy bands of FLGs become the dispersionless LLs. Some magneto-electronic properties 
of FLGs have been obtained by measurements of the quantum Hall conductivity, 2,20-22 
scanning tunneling spectroscopy, 22-24 and cyclotron resonance. 25 ' 26 The effective- mass ap- 
proximation (EM) and the tight-binding model (TB), two theoretical methods based on 
different concepts in calculations, are usually utilized to study the magneto-optical prop- 
erties of layered-graphene systems. In MG, EM shows that the linear bands become LLs 
characterized by the special relation E% v oc \^n c ' v B , 27 ' 28 where c (v) is the index of the 
conduction (valence) bands. E% v and n c ' v are the energy and the integer quantum number 
of the nth energy band respectively. However, TB indicates that the special relation will be 
broken as the energy exceeds the region of ±0.5 eV since the linear dispersions become grad- 
ually parabolic as energy exceeds this region. 29 As for the optical-absorption spectra, the 
low-frequency absorption peaks possess the selection rule |An| = |n c — n v \ = 1 and their 
frequency is proportional to \^B^, 30 which has been confirmed by magneto-transmission 
measurements. 25 ' 31 ' 32 
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For BBG, EM takes the LL structures of MG and uses them in the BBG calculations. 
Interlayer interactions are regarded as perturbations (EM considers only one or two in- 
terlayer interactions). Based on EM calculations, the first theoretical study 13 of LLs in 
BBG derived an effective two-dimensional Hamiltonian that describes the low-energy elec- 
tronic excitations. Furthermore, EM predicts the existence of two groups of the occupied 
and unoccupied LLs, and they are symmetric about Ep because the parabolic bands at 
Bq=0 are symmetric in the calculations. 30 ' 33 ' 34 The low-energy LLs in their results follow 
the relation E% v oc B . 30 As for the optical absorption spectra, EM calculations show 
that there are three categories of absorption peaks. 30 The absorption frequencies obey a 
linear function of B in weak fields, and depend linearly on v^o as the corresponding 
energy leaves the parabolic band region. 30 The three categories own the same optical se- 
lection rule | An | = i. 14 > 30 > 35 > 36 i n TB, a magnetic field and four most important interlayer 
interactions 8 ' 12 ' 13 are simultaneously included in the calculations. The magnetic flux 
through a hexagon corresponding to the strongest field strength (~60 T) in experiments 
is much smaller than the flux quantum O (0 ~ </> /1320). From the viewpoint of the TB 
perturbation theory, electronic states that lie in close proximity to each other at -Bo=0 
will be aggregated and become dispersionless LLs under the influence of a magnetic field. 
The characteristics of electronic states at B =0 are directly reflected in the main features 
of LLs. The interlayer interactions strongly affect the electronic properties, especially for 
low-energy levels. These interactions are neither negligible nor suitable as perturbations. 
TB indicates that the low-energy LLs should exhibit more complex features; 12 i.e., there 
exists an energy gap, two groups of asymmetric LLs ( ls *LL c,u s and 2nd LL c ' v s), and no simple 
relations between E% v and B . The main features would be reflected in the magneto-optical 
properties. 

In the initial stage of our works in studying the magneto-electronic properties, only 
eigenvalues and eigenvectors at strong magnetic field strength were obtained 37 because the 
Hamiltonian matrix is quite large, e.g., the matrix is about 31600x31600 for MG at 10 T. 
Through the band-like-matrix strategy 12,38-42 designed by our group, the eigenvalues and 
the wave functions at weak fields (~1 T) are thus clearly depicted. 12 However, the calcula- 
tion of optical properties is more complex and time-consuming than that of the electronic 
properties since detailed calculations for a huge velocity matrix and reliable characteriza- 
tion of the wave functions are necessary for the TB-based optical-absorption spectra. In 
this work, the localized features of the wave functions are utilized to effectively reduce 
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the computation time, and thus the optical absorption spectra can be solved. It's the 
first time that the optical properties of BBG calculated by TB are obtained. The influ- 
ence of the four most important interlayer interactions on the magneto-optical absorption 
spectra of BBG is investigated in detail by TB and gradient approximation. 37 > 4 °- 44 The 
wave functions are clearly depicted and used to divide the LLs into two groups, 12 lst LLs 
and 2nc LLs. The two groups lead to four categories of optical absorption peaks. These 
peaks exhibit twin-peak structures, complex selection rules, regular frequency-dependent 
absorption rates, and complicated field-dependent absorption frequencies. What deserves 
to be mentioned is that our study is more complete in nature and the presented results 
are more reliable than those of previous works. Previously experimental and theoretical 
research mainly focused on the optical properties associated with lst LLs. However, in our 
work, the optical absorption peaks resulting from both groups are discussed in detail and 
each peak can be clearly identified. A detailed comparison of optical-absorption spectra 
between MG and BBG is made. In addition, comparisons of optical properties between 
AA-stacked bilayer graphene (AABG) and BBG and between EM and TB are also made. 

RESULTS AND DISCUSSION 

The dispersionless LLs of BBG at B = 40 T are shown in Figure la. The interlayer 
interactions strongly influence the main features of LLs. The conduction and valence bands 
are asymmetric. There is an energy gap between the lowest unoccupied LL and the highest 
occupied LL. Based on the characteristics of wave functions, LLs can be classified into two 
groups, lst LLs (black curves) and 2nd LLs (red curves). ls *LLs appears at \E C ' V \ > 0, and 
2nd LLs begins at E c > 71 +7 6 (0.34 eV) and E v < -71 +7 6 (-0.39 eV). The onset energies 
of the first and second groups are consistent with those of the first and second pairs in the 
zero-field subbands. Furthermore, no simple relation between the energy and field strength 
can be found for these LLs, i.e., E% v is not directly proportional to either Bq or y/B~o. 12 

The characteristics of the wave functions deserve a closer investigation because they 
strongly dominate the optical properties. Each LL is composed of fourfold degenerate 
states with similar characteristics. At k x = k y = 0, the wave functions mainly consist of 
four localized subenvelope functions centered at positions X\ = 1/6, x 2 = 1/3, x 3 = 2/3, 
and rr 4 = 5/6. Xj = i_ 4 = A ijN /2R B or = B ijN /2R B . A i>N (B i>N ) indicates the Nth. A 
(B) atom on the ith layer, where N = 1,2...2R B and i — 1,2 denote the atomic site in 
the primitive unit cell and the layer index respectively (please see the detailed discussion 
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regarding these representations in the METHODS section). The feature of a wave function 
resulting from the atoms with an even number (Af* and B^) is similar to that with an 
odd (Afl and B°£) number, where e (o) is an even (odd) integer. Because of this, only- 
odd atoms are discussed in the following calculations. For convenience, one of the fourfold 
degenerate states at Xi = 1/6 is selected to be examined, as shown in Figures lb-le. 
Through appropriate fitting, the wave functions of the nith ls< LL (black curves) and n 2 th 
2nd LL (red curves) can be expressed as 

A c lo , A 2o oc ip 2 (x) , B[ oc ifi (x) , B c 2o oc (p (x) , n c { eff = for m = 1; (la) 
Al Q , A\ Q oc y? (x) , Bl Q oc y? 2 (x) , B\ oc Vl (x) , n\' eff = 1 for n x = 1; (lb) 
A%, A c 2 % oc (/V-i (x) , BZ oc <p ni _ 2 (x) , BZ oc v^ ni (x) for <' e// = m > 2. (lc) 



- A i',o> ^2,0 oc (x) , oc (p 2 (x) , 5^ oc ^ (x) , n^' 6 ^ = for n 2 = 1; (2a) 
^f , ^f oc ^n 2 -i (x) , SJJ oc ^„ 2 - 2 (x) , 5 2 c ;: oc v^„ 2 (x) for nr' e// = n 2 - 1 > l(2b) 

(x) is the product of the nth-order Hermite polynomial and Gaussian function, 12 ' 29 ' 40 
where n is the number of zero points of ip n (x) and chosen to define the quantum number 
of a LL. 12 That is to say, there are four quantum numbers in each LL since a LL wave 
function is made up of four magnetic TB functions. Thus an effective quantum number, 
n e ff, defined by one of these four quantum numbers is necessary to identify each LL. 
The effective quantum number of a LL in each group is assigned by the number of zero 
points in one of the four magnetic TB functions, with the largest amplitude among these 
four TB functions at the onset energy of each group. Accordingly, n c { v,e (n^ ) is the 
effective quantum number of the ls *LLs ( 2nd LLs) defined by the quantum number of B 2 V 
( A^o)- It should be noted that the wave functions of two groups are composed of not 
only the B atoms but also the A atoms, which are different from those calculated by EM 
including only contributions from the B atoms. 13 ' 14 Such difference could cause the optical 
properties obtained by TB and EM to exhibit different features. lst LLs and 2rad LLs can be 
distinguished by the different features of the wave functions in the mixed region (E c > 0.34 
eV and E v < —0.39 eV), which is useful in analyzing the optical excitation channels. 

When BBG is subjected to an electromagnetic field at zero temperature, only excitations 
from the occupied to the unoccupied bands occur. Based on Fermi's golden rule, the optical- 
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absorption function is given by 

a(u) oc Yl I 7^\mnifM—\nnii ff M\ 

r eff „ Rtt JistBZ W rn e 
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c,eff v,eff 
n l,2 i n l,2 



xlm{ L , V N Ji , V \ ^}, (3) 

(k, n-f) - E v (k, nlf f ) - u - iT 

where E is the unit vector of an electric polarization. T (=1 meV) is a broadening parameter 
and often affected by temperature and defect effects. If a purer sample is manufactured and 
observed its physical properties under a sufficiently low temperature, T will be small enough 
for observing important fine structures. 22 ' 45 Furthermore, (ip c (nl 2 , k)|E ■ ~P / m e \ijj v (n\ 2 , k)) 
is the velocity matrix element derived from the dipole transition and denoted M cv . In this 
work, M cv is calculated by gradient approximation. Through detailed calculations, M cv 
can be simplified to be the product of three matrices, the initial state (occupied state; 
t/j v ), the final state (unoccupied state; tp c ) and VkHij. The last term corresponds to the 
direction of electric polarization, i.e., it is dHij/dk x for polarization along the x-axis in 
this work. The elements of dHij/dk x are non-zero only for associated with the hopping 
integrals, as are the velocity matrix elements. In other words, when the velocity matrix 
element does not vanish, the initial and final states in the product should be the two states 
corresponding to the non-vanishing hopping integrals. Accordingly, M cv is expressed as 

2R B 

{ E E ^ A io + A 'eY >< i B l' + *&)]V fc (A,JVk Wb\ B hN , k ) 
i=l,2 N,N'=1 

2R B 

+ E E + x ( B ]o> + B] e ,)]V k (B hNk \H B \ B jtNnt ) 

i,j=l,2;ijij N,N'=1 
2R B 

+ E E K S * + BIT * {A* + A )e)Wk (B hNk \H B \ A jiN , k )} + h.c. (4) 

i,j=l,2;i=ij N,N'=1 

The first, second, and third term in eq 4 correspond to the three hopping integrals 70, 
73, and 74, respectively. That is to say, the optical excitations are provided by the 
three channels related to 70, 73 and 74. For convenience, (A c io + A c ie )* x [B^ , + ££,,), 
(S? + BIT x (B] , + B] e ,), and (S? + Bf e )* x + Aj e ,) are represented by M^( 7o ), 
^^3(73), and (74), respectively in the following discussions. It should be noted that 
eq 4 implies that the characteristics of the wave functions would play a dominant role in 
determining the selection rule and absorption rate in optical excitations. 

The low-frequency optical-absorption spectra of BBG regarding four different magnetic 
field strength are shown in Figures 2a-2d, respectively. The absorption spectrum of MG 
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is also shown in Figure 2e. In Figure 2a, the absorption spectrum of BBG at B = 
40 T presents prominent and inconspicuous (indicated by the black arrows) peaks. The 
prominent peaks with definite optical selection rules are discussed in detail in the following 
paragraphs. The peaks can mainly be classified into four categories of peaks, Wn's (black 
dots), ou 22 s ( r ed dots), ou^s (green dots), and u; 2 i's (blue dots), which originate in the four 
different excitation channels, lst LL v s to lst LL c s, 2nd LL v s to 2nd LL c s, lst LL"s to 2nd LL c s, and 
2nd LL v s to ls *LL c s respectively. The former (latter) two channels result from the transitions 
between two LLs in the same (different) groups, as illustrated by the inset of Figure 2d, 
and they display obvious optical selection rules. 

The optical excitations associated with each prominent peak can clearly be identified, 
as shown in Figure 2a. For convenience, the excitations of Wn's, ou 22 s, ou i2 s, and w 2 i's are 
respectively represented as n{ e ^ — > n{ e ^ , n 2 e ^ — > n 2 e ^ , n v { e ^ — > n 2 e ^ , and n 2 e ^ — > 
n ^ e ff j n ^he f n ow i n g. The selection rule of cu^'s is denoted by An^ (= n c - e ^ — ). As 
to Wn's and W22's, only the first peak (cuLj of Wn's originates in single channel 1 — > 2. The 
other mth peak (w^) consists of the pair and ou™' H (uj™ 2 L and io™ 2 ) corresponding 
to m — > m + 1 and m + 1 — > m (m— 1 — > m and m — >■ m — 1) respectively. The superscript 
L (H) indicates the lower (higher) energy of the pair. For example, u; n (^22) * s composed 
of w n L and u 2 ^ {u\ 2 and <jo 22 h ) owing to 2 — > 3 and 3 — > 2 (0 -> 1 and 1 — )■ 0) respectively. 
The origin of a twin-peak structure is that the energies of m — > m + 1 (m — 1 — > m) and 
m + 1 — )• m (m — > m — 1) in Wn's (w 22 's) are slightly different owing to the asymmetry of 
the occupied and unoccupied LLs. It should be noted that similar twin-peak structures are 
also shown in a recent experiment conducted to study the magneto-absorption spectra. 45 In 
short, the optical selection rules of Wn's and oj 22 s can be represented by Ann = An 22 = 
±1, which is the same as the LLs in MG. 

As for the intergroup transitions, twin-peak structures are also obtained in ui\ 2 s and 
w 2 i's. The mth peak of Wi 2 's is formed with the pair, u™ 2 L originating in m — > m and 
lo™ 2 h resulting from m + 1 — > m — 1. For instance, the first pair is indicated by io\ 2 and 

1 H 7TL L 

oo 12 in Figure 2a, which correspond to 1 — > 1 and 2 — > respectively. That is to say, oo 12 
owns the optical selection rule Ani 2 = and uf^f 1 possesses An 12 = —2. For w 2 i's, the 
excitations — > and — > 2 lead to the peaks u 21 and u 2l of the first pair in Figure 2a 
respectively. The other channels m — 1 — > m + 1 and m — > m result in the pair u^' 1 and 
0J 2i H respectively. In other words, the selection rules of w 2 i's are An 2 i = and An2i = 2. 
Obviously, the optical selection rules of Ui 2 s and u; 2 i's are different from those of Wn's and 
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cu 22 's, a fact which is not clearly discussed in the previous theoretical works. 

The optical selection rules can be explained from two different aspects, i.e., the effective 
quantum numbers and the numbers of zero points between the initial and final states. 
Both aspects are based on the premise that the velocity matrix is derived from the dipole 
transition. The product of the conduction and valence wave functions in eq 4 would not 
only determine the allowed excitation channel, but also enable us to estimate the absorption 
rate. This means that the dissimilarities in selection rules of the four categories mainly 
result from the characteristics of the wave functions. In eq 4, M^ b {^q) associated with 7 
dominates the excitations of the prominent peaks. It strongly depends on the numbers of 
zero points of A c { v and B c { v (A% v and B c 2 ' v ) since MJ B ( 7o ) = {A c io + A c ie )* x {B v io , + £?,)■ 
-^ab(To) nas non-zero values only when A c { v and B c { v {A c 2 ,v and B 2 V ) own the same <£ n (x) 
because of the orthogonality of (f n (x). Thus, the selection rules of prominent peaks can be 
easily obtained by the definition of the wave functions in eqs I and 2. From the aspect of 
the numbers of zero points between Af v and B° ,v , these states own the same number of zero 
points in an available excitation channel in both intragroup and intergroup transitions. As 
for the effective quantum numbers, the selection rules are Ann = An 22 = ±1, An 12 = 0, 
—2, and An 2 i = 0, 2. Although these selection rules are different, all excitations obey the 
underlying condition that the numbers of zero points of the initial and final states in the 
dipole-transition-derived velocity matrix are the same. 

In addition to the optical selection rules, the peak intensity and absorption frequency 
are also discussed in this work. The absorption rate significantly relies on the field strength 
and the excitation energy. In Figures 2a-2d, the peak intensity is increased with enlarging 
field strength, while the opposite is true for the peak number. The peak intensity under 
a fixed magnetic field strength is determined by the product of two wave functions with 
respect to the two nearest-neighbor atoms. The wave functions of ls *LLs ( 2nd LLs) are 
mainly dominated by the B (A) atoms. The amplitudes of A /" and B°' v , respectively, 
gradually grow and decrease (decrease and grow) with increasing energy in lsi LLs ( 2nd LLs). 
Accordingly, the value of the product of Af v and B°' v increases (decreases) with rising 
energy for w u 's and w 22 's (wi 2 's and w 2 i's). This value is largest when the amplitudes of Af v 
and B° ,v are equal. The intensities of w n 's and w 22 's (o>i 2 's and w 2 i's) thus increase (decline) 
with increasing frequency. Furthermore, some peaks with commensurate energies overlap 
because the peak spacing declines as the frequency increases; this overlap is responsible for 
the variation in peak intensity. In contrast to BBG, MG exhibits absorption peaks with 
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uniform intensity (Figure 2e). 

The field-dependent absorption frequencies (a> a 's) of u u 's, Co'22's, ou^s and W2i's are 
shown in Figure 3a by the black, red, green, and blue dots, respectively. w a 's of the four 
categories rise with increasing B . These absorption frequencies do not exhibit the simple 
relations indicated by EM calculations, 29 i.e., the frequencies are not directly proportional 
to either B or v^o- Apparently, oo a oc v^o hi MG (Figure 3b) is very different from 
that in BBG. The main reason for the difference between BBG and MG is that interlayer 
interactions play important roles in BBG and are not negligible in calculations. In addition, 
the convergent frequencies of Wn's, c^'s, Wi 2 's and w 2 i's at the weak field strength are 
approximately 0, 0.73 eV (271), 0.34 eV (71 + 76), and 0.39 eV (71 — 76) respectively. This 
implies that the optical measurements can reasonably determine the values of 71 and 7 6 
through observing the convergent frequencies of absorption peaks. The predicted results 
are very useful and reliable for experimental studies. 

The inconspicuous peaks indicated by arrows in Figure 2a also obey specific selection 
rules (not shown). These peaks mainly originate from excitations related to two interlayer 
interactions (73 and 74) and differ from the prominent peaks associated with the intralayer 
hopping integral (70)- The inconspicuous peak intensities are much weaker than the inten- 
sities of the prominent peaks since the values of the velocity matrix elements corresponding 
to 73 and 74 are much smaller than those related to 70 . The different origins and selection 
rules of the prominent and inconspicuous peaks imply that the optical properties can not 
be simply derived from the eigenenergies of LLs. The detailed calculations of the velocity 
matrix and reliable characterization of the wave functions are necessary to determine the 
optical selection rules and the peak intensity. 

The atomic hopping integrals are important to analyze the physical properties of BBG. 
So far the values of these atomic interactions have not been well defined. Figure 4a is the 
optical-absorption spectrum calculated with the parameters used in our work. 8 Figures 4b, 
4c, and 4d are the results calculated with three other sets of parameters 10,16,17 used in the 
investigations of BBG, respectively. The optical-absorption spectra corresponding to the 
four sets of parameters show similar qualitative properties, i.e., they exhibit four categories 
of absorption peaks, twin-peak structures, and same optical selection rules. However, the 
absorption frequencies in the four results are slightly different since the values of the four 
sets of parameters are distinct. The reliability of these parameters could be verified through 
more accurate experimental measurements with a purer bilayer graphene. 
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TB and EM methods are based on different concepts in the BBG calculations. The 
former simultaneously includes a magnetic field and four important interlayer interactions 
in the calculations. However, the latter takes the LL structures of MG to study the phys- 
ical properties of BBG, where interlayer interactions are regarded as perturbations. The 
wave functions in TB calculations are well characterized and different from those in EM 
calculations. The magneto-optical absorption spectra studied by TB and EM exhibit five 
crucial differences. The former demonstrate four categories of absorption peaks, twin-peak 
structures in each category, complex optical selection rules, and complex field-dependent 
absorption frequencies. However, the latter display three peak categories, no twin-peak 
structure, a single selection rule, and the simple field-dependent absorption frequencies 
uj a oc B and u a oc ^/Bq. Moreover, the regular frequency-dependent peak intensities are 
shown in TB results, while they are not obtained in EM results. These differences mainly 
originate from the different concepts and distinct definitions of the wave functions in the 
two methods. 

In addition to BBG, AABG also presents special optical properties, as shown in Figures 
S6 and S7 of the supporting material. AABG and BBG show similar LL features, i.e., two 
groups of LLs and an asymmetric structure. However, these two bilayer graphenes display 
different optical properties. AABG contains two categories of absorption peaks resulting 
from the intragroup excitations and these peaks exhibit only one single optical selection 
rule. No twin-peak structure similar to that in the BBG case arises from the special relations 
in the two LL groups of AABG. The differences of the optical absorption spectra of AABG 
and BBG imply that the optical properties can not be directly identified through only the 
eigenvalues of LLs. An explicit analysis of the characteristics of the wave functions and 
detailed calculations of the velocity matrix elements are necessary to calculate the optical 
properties. Both the magneto-electronic and magneto-optical properties of AABG and 
BBG reflect the main features of zero-field electronic and optical properties. A detailed 
analysis of the magneto-optical absorption spectra in AABG is discussed in an unpublished 
work (Y. H. Ho et al, Characterization of Landau-Level Wavefunctions and Magneto- 
Optical Excitations in AA-stacked Bilayer Graphene). 

TB is widely applied to tackle physical problems of semiconductors and carbon-related 
systems. Through the strategy developed by our group, the complex Hamiltonian matrix 
can be solved. The eigenvalues and wave functions are clearly depicted in this calculation, 
and thus can be used to interpret the results in optical properties. Since TB calculations 
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simultaneously include the magnetic field and the important interlayer interactions, such 
results are more reliable and accurate over a much wider frequency range. Besides, TB 
can conceivably be expanded to investigate other layered-graphene systems (multilayer 
graphenes and bulk graphite) and other physical properties (electronic excitations and 
transport properties). 

SUMMARY AND CONCLUSIONS 

In summary, it can be said that the interlayer interactions play very important roles in 
the magneto-optical properties of BBG. In this work, detailed calculations of the velocity 
matrix and reliable characterization of the wave functions are provided for analyzing the 
optical properties. Based on the characteristics of the wave functions, the asymmetric LLs 
can be divided into two groups. These two groups lead to four categories of prominent 
optical-absorption peaks Wn's, W22'S) Wi 2 's, and w 2 i's. Each category exhibits twin-peak 
structures originating from the asymmetry of LLs. The optical selection rules of the four 
categories are Ann = ±1, Ari22 = ±1, Ani2 = 0,-2, and An 2 i = 0,2 respectively. All 
of the selection rules can be reasonably explained by the characteristics of the wave func- 
tions. The peak intensities of Wu's and w 2 2's (^12's and w 2 i's) increase (decrease) with 
rising frequency. The field-dependent absorption frequencies are not directly proportional 
to either B or ^/Bq. Furthermore, the convergent absorption frequencies at the weak 
field strength might be helpful and reliable in determining the interlayer atomic interac- 
tions 71 and 76. The main results of BBG in this work are more complex in nature than 
those of MG, and different from those in previously published works studied by EM. The 
above-mentioned magneto-optical properties could be confirmed by magneto-absorption 
spectroscopy measurements. 26,31 ' 32 

METHODS 

The optical-absorption spectrum is calculated by gradient approximation based on the 
tight-binding model. The geometric configuration of a primitive unit cell in BBG is shown in 
Figure 5a. The unit cell includes four sublattices, A ± , B 1 , A 2 and B 2 . The subscripts 1 and 
2 are, respectively, the indices of the first and second layer. The nearest-neighbor hopping 
integral on the same layer is indicated by 7 (=2.598 eV). The four important interlayer 
interactions 8 used in this work are 71 (=0.364 eV), 73 (=0.319 eV), 74 (=0.177 eV), and 7 6 
(=-0.026 eV), indicated in Figure 5a. 7 6 is the difference of site energy between the A and B 
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atoms owing to the interlayer coupling. The perpendicular-uniform magnetic field induces 
a periodic Peierls phase related to the vector potential A = (0, B x, 0). The primitive unit 
cell shown in Figure 5b is enlarged under this periodic condition and owns four effective 
sublattices including 2R B Ai, 2R B Bi, 2R B A 2 , and 2R B B x atoms, respectively. R B is 
related to the dimension of the Hamiltonian under a magnetic field; 12 for example, R B is 
7900 for Bo = 10 T. That is to say, each LL is the linear combination of the four magnetic 
TB functions associated with the four effective sublattices. The representation A^x (B^n) 
indicates the Nth (N = 1, 2...2R B ) A (B) atom on the ith (i = 1, 2) layer. The Hamiltonian 
matrix elements in the presence of a magnetic field are 

(R-i,Af|-^B|Rj',Af') = ls(Ri,N, R-i',AT') ^ exp[?k • (Ri',jV' — Rj,jv) + Z— AG^R^, Ri^JV 7 )]) 

(5) 

where R^jv is the position vector of A itN or B i N . 7 s 's(R i> Ar, Rj' t jv) indicate the atomic 
interactions between the atoms at R^jv an d Ri',jv'> *- e -> they are 70, 71, 73, 74, and j e . 
AG(Rj ) Ar, Rj',jv') = fo(Ri',Ni — Ri.at) • AfRj/^/ + A(Ri' i jv — TU,N)]dX is the Peierls phase 
induced by a magnetic field (Please see more detailed definitions of the Hamiltonian matrix 
elements in Ref. 12). After diagonalizing the Hamiltonian matrix, the eigenvalues are 
obtained and thus the optical-absorption function in eq 3 can be solved through gradient 
approximation. 

Acknowledgments. This work was supported by the NSC and NCTS of Taiwan, under the 
grant Nos. NSC 95-2112-M-006-028-MY3 and NSC 97-21 12-M-110-001-MY2. 

Supporting Information Available: Influences of the interlayer interactions on the magneto- 
electronic and magneto-optical properties (Figures SI and S2), Landau level and magneto- 
optical absorption spectra calculated with four sets of parameters (Figures S3 and S4), opti- 
cal excitations obtained with four broadening parameters (Figure S5), and electronic struc- 
tures and optical-absorption spectra (Figures S6 and S7) of AA-stacked bilayer graphene. 



This material is available free of charge via the Internet at http://pubs.acs.org 



12 



References 



1 Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; 
Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 
306, 666. 

2 Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Katsnelson, M. L; Grig- 
orieva, I. V.; Dubonos, S. V.; Firsov, A. A. Two-Dimensional Gas of Massless Dirac 
Fermions in Graphene. Nature 2005, 438, 197-200. 

3 Berger, C; Song, Z. M.; Li, T. B.; Li, X. B.; Ogbazghi, A. Y.; Feng, R.; Dai, 
Z. T.; Marchenkov, A. N.; Conrad, E. H.; First, P. N.; de Heer, W. A. Ultrathin 
Epitaxial Graphite: 2D Electron Gas Properties and a Route toward Graphene-Based 
Nanoelectronics. J. Phys. Chem. B 2004, 108, 19912-19916. 

4 Rutter, G. M.; Crain, J. N.; Guisinger, N. P.; Li, T.; First, P. N.; Stroscio, J. A. 
Scattering and Interference in Epitaxial Graphene. Science 2007, 317, 219-222. 

5 Wallace, P. R. The Band Theory of Graphite. Phys. Rev. 1947, 71, 622-634. 

6 Griineis, A.; Attaccalite, C; Wirtz, L.; Shiozawa, FL; Saito, R.; Pichler, T.; Rubio, A. 
Tight-Binding Description of the Quasiparticle Dispersion of Graphite and Few-layer 
Graphene. Phys. Rev. B 2008, 78, No. 205425. 

7 Slonczewski, J. C; Weiss, P. R. Band Structure of Graphite. Phys. Rev. 1958, 109, 
272-279. 

8 Charlier, J.-C; Michenaud, J. -P. First-Principles Study of the Electronic Properties 
of Simple Hexagonal Graphite. Phys. Rev. B 1992, 46, 4531-4539. 

9 Latil, S.; Henrard, L. Charge Carriers in Few-Layer Graphene Films. Phys. Rev. 
Lett. 2006, 97, No. 036803. 

Partoens, B.; Peeters, F. M. From Graphene to Graphite: Electronic Structure around 
the K Point. Phys. Rev. B 2006, 74, No. 075404. 

1 Lu, C. L.; Chang, C. P.; Huang, Y. C; Chen, R. B.; Lin, M. F. Influence of an 
Electric Field on the Optical Properties of Few-Layer Graphene with AB Stacking. 
Phys. Rev. B 20 06, 73, No. 144427. 



13 



Lai,Y. H.; Ho, J. H.; Chang, C. P.; Lin, M. F. Magnetoelectronic Properties of Bilayer 
Bernal Graphene. Phys. Rev. B 2008, 77, No. 085426. 

McCann, E.; Fal'ko, V. I. Landau- Level Degeneracy and Quantum Hall Effect in a 
Graphite Bilayer. Phys. Rev. Lett. 2006, 96, No. 086805. 

Abergel, D. S. L.; Fal'ko, V. I. Optical and Magneto-Optical Far- Infrared Properties 
of Bilayer Graphene. Phys. Rev. B 2007, 75, No. 155430. 

Charlier, J.-C; Gonze, X.; Michenaud, J. -P. First- Principles Study of the Electronic 
Properties of Graphite. Phys. Rev. B 1991, 43, 4579-4589. 

Nakao, J. Landau Level Structure and Magnetic Breakthrough in Graphite. J. Phys. 
Soc. Japan 1976, 40, 761-768. 

Kuzmenko, A. B.; Crassee, I.; van der Marel, D. Determination of the Gate- Tunable 
Band Gap and Tight-Binding Parameters in Bilayer Graphene Using Infrared Spec- 
troscopy. Phys. Rev. B 2009, 80, No. 165406. 

Ohta, T.; Bostwick, A.; Seller, T.; Horn, K.; Rotenberg, E. Controlling the Electronic 
Structure of Bilayer Graphene. Science 2006, 313, 951-954. 

Ohta, T.; Bostwick, A.; McChesney, J. L.; Seyller, T.; Horn, K.; Rotenberg, E. 
Interlayer Interaction and Electronic Screening in Multilayer Graphene Investigated 
with Angle- Resolved Photoemission Spectroscopy. Phys. Rev. Lett. 2007, 98, No. 
206802. 

Novoselov, K. S.; McCann, E.; Morozov, S. V.; Fal'ko, V. L; Katsnelson, M. L; 
Zeitler, U.; Jiang, D.; Schedin, F.; Geim, A. K. Unconventional Quantum Hall Effect 
and Berry's Phase of 2vr in Bilayer Graphene. Nat. Phys. 2006, 2, 177-180. 

Zhang, Y.; Tan, Y.-W.; Stormer, H. L.; Kim, P. Experimental Observation of the 
Quantum Hall Effect and Berry's Phase in Graphene. Nature 2005, 438, 201-204. 

Miller, D. L.; Kubista, K. D.; Rutter, G. M.; Ruan, M.; de Heer, W. A.; First, P. 
N.; Stroscio, J. A. Observing the Quantization of Zero Mass Carriers in Graphene. 
Science 2009, 324, 924-927. 

Li, G.; Andrei, E. Y. Observation of Landau Levels of Dirac Fermions in Graphite. 
Nat. Phys. 2007, 3, 623-627. 

14 



Matsui, T.; Kambara, H.; Niimi, Y.; Tagami, K.; Tsukada, M.; Fukuyama, H. STS 
Observations of Landau Levels at Graphite Surfaces. Phys. Rev. Lett. 2005, 94, No. 
226403. 

Deacon, R. S.; Chuang, K.-C; Nicholas, R. J.; Novoselov, K. S.; Geim, A. K. Cy- 
clotron Resonance Study of the Electron and Hole Velocity in Graphene Monolayers. 
Phys. Rev. B 20 07, 76, No. 081406. 

Henriksen, E. A.; Jiang, Z.; Tung, L.-C; Schwartz, M. E.; Takita, M.; Wang, Y.-J.; 
Kim, P.; Stormer, H. L. Cyclotron Resonance in Bilayer Graphene. Phys. Rev. Lett. 
2008, 100, No. 087403. 

McClure, J. W. Diamagnetism of Graphite. Phys. Rev. 1956, 104, 666-671. 

Zheng, Y.; Ando, T. Hall Conductivity of a Two-Dimensional Graphite System. Phys. 
Rev. B 2002, 65, No. 245420. 

Ho, J. H.; Lai, Y. H.; Chiu, Y. H.; Lin, M. F. Landau Levels in Graphene. Physica 
E 2008, 40, 1722-1725. 

Koshino, M.; Ando, T. Magneto-Optical Properties of Multilayer Graphene. Phys. 
Rev. B 2008, 77, No. 115313. 

Sadowski, M. L.; Martinez, G.; Potemski, M.; Berger, C; de Heer, W. A. Landau 
Level Spectroscopy of Ultrathin Graphite Layers. Phys. Rev. Lett. 2006, 97, No. 
266405. 

Plochocka, P.; Faugeras, C; Orlita, M.; Sadowski, M. L.; Martinez, G.; Potemski, 
M.; Goerbig, M. O.; Fuchs, J.-N.; Berger, C; de Heer, W. A. High-Energy Limit of 
Massless Dirac Fermions in Multilayer Graphene Using Magneto- Optical Transmis- 
sion Spectroscopy. Phys. Rev. Lett. 2008, 100, No. 087401. 

Benfatto, L.; Sharapov, S. G.; Carbotte, J. P. Robustness of the Optical Conductivity 
Sum Rule in Bilayer Graphene. Phys. Rev. B 2008, 77, No. 125422. 

Nicol, E. J.; Carbotte, J. P. Optical Conductivity of Bilayer Graphene With and 
Without an Asymmetry Gap. Phys. Rev. B 2008, 77, No. 155409. 



15 



Mucha-Kruczyhski, M.; Abergel, D. S. L.; McCann, E.; Fal'ko, V. I. On Spectral 
Properties of Bilayer Graphene: The Effect of an SiC Substrate and Infrared Magneto- 
Spectroscopy. J. Phys.: Condens. Matter 2009, 21, No. 344206. 

Mucha-Kruczyrlski, M.; McCann, E.; Fal'ko, V. I. The Influence of Interlayer Asym- 
metry on The Magnetospectroscopy of Bilayer Graphene. Solid State Communica- 
tions 2009, 149, 1111-1116. 

Chang, C. P.; Lu, C. L.; Shyu, F. L.; Chen, R. B.; Fang, Y. K.; Lin,M. F. Magneto- 
electronic properties of a graphite sheet. Carbon 2004, 42, 2975-2980. 

Chiu, Y. H.; Lai, Y. H.; Ho, J. H.; Chuu, D. S.; Lin, M. F. Electronic Structure of 
a Two- Dimensional Graphene Monolayer in a Spatially Modulated Magnetic Field: 
Peierls Tight-Binding Model. Phys. Rev. B 2008, 77, No. 045407. 

Ho, J. H.; Lai, Y. H.; Chiu, Y. H.; Lin, M. F. Modulation Effects on Landau Levels 
in a Monolayer Graphene. Nanotechnology 2008, 19, No. 035712. 

Huang, Y. C; Chang, C. P.; Lin, M. F. Magnetic and Quantum Confinement Effects 
on Electronic and Optical Properties of Graphene Ribbons. Nanotechnology 2007, 
18, No. 495401. 

Chiu, Y. H.; Ho, J. H.; Chang, C. P.; Chuu, D. S.; Lin, M. F. Low-Frequency 
Magneto-Optical Excitations of a Graphene Monolayer: Peierls Tight-Binding Model 
and Gradient Approximation Calculation. Phys. Rev. B 20 08, 78, No. 245411. 

Huang, Y. C; Lin, M. F.; Chang, C. P. Landau Levels and Magneto-Optical Prop- 
erties of Graphene Ribbons. J. Appl. Phys. 2008, 103, No. 073709. 

Lin, M. F.; Shung, K. W.-K. Plasmons and Optical Properties of Carbon Nanotubes. 
Phys. Rev. B 1994, 50, 17744-17747. 

Johnson, J. G.; Dresselhaus, G. Optical Properties of Graphite. Phys. Rev. B 1979, 
7, 2275-2285. 

Chuang, K.-C; Baker, A. M. R.; Nicholas, R. J. Magnetoabsorption Study of Landau 
Levels in Graphite. Phys. Rev. B 2009, 80, No. 161410(R). 



16 



FIGURE CAPTIONS 



Figure 1. (a) Landau levels of bilayer Bernal graphene at B =A0 T. The wave functions 
of (b) Ai t0 , (c) Bi t0 , (d) A 2t0 , (e) -82,0 atoms with odd integer indices are shown. 

Figure 2. Optical absorption spectra of bilayer Bernal graphene at (a) 40 T, (b) 30 T, 
(c) 20 T, and (d) 10 T. The spectrum of monolayer graphene at 40 T is plotted in (e). 

Figure 3. Field-dependent optical absorption frequencies of (a) bilayer Bernal graphene 
and (b) of monolayer graphene. 

Figure 4. (a) The optical-absorption spectra at B = 40 T calculated with the param- 
eters used in our work, (b), (c), and (d) are the results calculated with three other sets of 
parameters commonly used in the investigations of bilayer Bernal graphene. 

Figure 5. The primitive unit cells of bilayer Bernal graphene in the absence and presence 
of a magnetic field are shown in (a) and (b) respectively. 70 (=2.598 eV) is the nearest- 
neighbor hopping integral and the four important interlayer interactions are 71 (=0.364 
eV), 73 (=0.319 eV), 74 (=0.177 eV), and 7 6 (=-0.026 eV). 
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